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ABSTRACT ~~ 


The "perpetual salt fountain" can occur in an area of the ocean 
where the water is vertically stable, Warm, saline water overlies cold, 
relatively fresh water. A tube is placed in the water from just beneath 
the ocean surface to just beneath the halocline. Water is pumped up the 
tube to initiate the flow. As the water in the tube is warmed, it be- 
comes less dense than (and thus buoyant with respect to) the surrounding 
water, and the flow will continue after pumping is stopped. 

The effect of the vertical conduction of heat by turbulent mixing 
on the velocity of flow in the "perpetual salt fountain" is investigated. 
Temperature and salinity profiles with depth from hydrographic soundings 
are considered to determine flow velocities in the ocean, 

The calculations show that vertical heat conduction has a negligible 
effect on the rate of flow in the presence of horizontal heat conduction, 
the value of the average salinity of a column of water minus the salinity 
of the water entering the bottom of the tube is the primary factor in- 
fluencing the rate of flow, and the effect of the temperature-depth profile 
on the velocity of flow is negligible so long as initial stability re- 


quirements are met. 
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1. Introduction. 

Stommel, Arons, and Blanchard [8] first wrote about an "oceano- 
graphic curiosity" which they called the "perpetual salt fountain." 
Consider a situation in which there is a layer of cold, relatively 
fresh water which underlies a wam, saline upper layer of water. The 
density of water is a function of temperature, salinity and pressure, 
and the density of the upper layer must be less than the density of 
the lower layer for stability. 

Now insert a tube into the water which extends from just beneath 
the surface of the water to the top of the lower layer of water. A 
power source is provided to initiate an upward flow of water. As the 
water in the lower layer is started up the tube, it is warmed by heat 
flowing into the tube from the surrounding warm water. Since the 
water in the tube is less saline than the water outside the tube, it 
becames less dense than the outside water and is forced upward. 

The flow process requires no source of power other than that used 
to initiate the flow, and flow will continue so long as there continues 
to be a well-defined two-layer system. The principle will work for a 
flow in either direction, 

Groves [2] calculated the rate of this flow in the tube to de- 
termine what, if any, practical value this perpetual flow phenomenon 
might have, Groves' calculation was based on a simple model where the 
water temperature was considered a constant down to the bottom of the 
tube which just extended into water of a colder temperature. Groves 
determined that the volume rate of flow for his tube would be too small 
to be of practical significance for fish farming (that is, use of this 


principle as a source for providing nutrient-rich deep water at the 


surface to increase the amount of edible fish in an area). 


Groves' model neglected the effects of vertical heat conduction 
within the tube, which should increase the warming of the water inside 
the tube and provide a corresponding increase in rate of flow. If this 
were the case, then perhaps a larger tube could be used. 

The additional heating which is a result of vertical heat con- 
duction would be sufficient to warm the increased volume of water in 
a larger tube. If the flow velocity in a tube of larger inner radius was 
approximately equal to the flow velocity in the tube of smaller inner 
radius, the tube with the larger inner radius would provide an increase 
in volume rate of flow. 

It was considered that a significant increase in the velocity of 
the water in the tube might increase the volume of flow sufficiently to 
warrant further investigation of this process and its practical appli- 
cations. The effects of wall friction on the temperature of the water 
inside the tube were considered, but the result was so small that the 
effect was considered negligible. Calculations of temperature change 
due to this stress are contained in Appendix C, 


2. Solution for a model which includes the effects of vertical 
eat conduction, 


In order to see the effect of vertical heat conduction on the 
rate of flow of the water inside the tube, let us consider the same 
model that Groves worked with, but include vertical heat conduction. 
Assume that the water density is a linear function of temperature 
and salinity, and that the water is incompressible. Therefore, the rate 
of change of density with temperature, de/ye , and the rate of change 
of density with salinity, °°/yS » are constant, a valid assumption 


within the small range of temperature and salinity considered. Then, 


e = Po + e (*%o), + §S (%s), (1) 


where Q is the water density and © and S respectively are the temp- 
erature and salinity of the water, The heat flow into the tube is 


calculated by the equation 


=f (e+? dy) =v 211+ dy 





=, ®, KT + dy a2 (2) 
where D/Dt designates the material derivative following the mean 
motion of the water, q is the heat per unit mass, n is the heat con- 
ductivity of tube material ( say, copper ), r is the inner radius of 
the tube, Z(y) is the temperature of the water in the tube, w is the 
wall thickness of the tube, and K is the coefficient of vertical heat 


conduction for turbulent flow calculated in Appendix D. At the sea 


surface y=0, and y is positive downward. If we assume a steady state, 
DO _ a0 3% __ 3 
Dt oy gt JY (3) 
since the local change an =O. ( The average flow velocity is v.) 
Cambining the -_" tems, equation (2) becomes 
do +g -«p=--xe (4) 
a and X=2n/Krw. The specific heat of sea water is c, h is 
the depth of the bottom of the tube, and L is the length of the tube. 


The general solution of the homogeneous equation associated 


with (4) is 


~, (y) = Tet + Jet (5) 


where J, and Jz are undetermined constant coefficients, 


i -6- ( 8" +420) 


2 (6) 
» _- B+ (Rs 4x) (7) 
2 oe & 


2.1. Solutions for fourth-power temperature and salinity profiles. 


In order to find a particular solution to equation (4), a fourth- 


power temperature profile was assumed: 


Oly) = A-Cy?, (8) 


This function was chosen because it can be easily fit to actual 
temperature-depth curves, The value of sea-surface temperature can 
be inserted»for A, andthe value (T;- T,,) /h" inserted for C. The 


complete solution for equation (4) is 


O(y)= i ed + J,e"# +A- cir [2ca? 129248] 


- 4 |i + 248%] 4 rs =, fizca+i2e%c/ y? 
age PR ey? - Cols 


The water entering the bottom of the tube is at the same temper- 


(9) 


ature as the surrounding water at that depth; therefore, §(h)#0(h). 

For steady flow, the water at the top of the tube is assumed to be at 
the same temperature as the water surrounding it. This is a reasonable 
assumption since we are considering vertical heat conduction, and the 
water at the top of the tube undergoes vertical turbulent mixing with 
the surrounding water. Thus, the other boundary condition is 
9(0)=0(0). Only upward flow is considered in this paper. 


Applying the boundary conditions to equation (9) and solving for 


Jandel 
(4) = [pore hepene oA 
Oly awh pazh 
p(erh-}) - Eh- Fh? -Gh3 o 424 (10) 
C inhle ah 


+ A- D-Ey-Fy*-Gy?-Cy? o 


The constants D, E, F, and G are functions of Xk, @ , and C for a given 


velocity: 
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d= ty [24cu? +728°CH + 248 ‘o], 
E= 25 /4Pecx + 24 8?c], 
Fe yliac% +12 ez c] ; 
G- alec]. 


Now assume that the water is removed fran the top of the tube as 


il 


fast as it flows up so that no additional pressure head is formed. 

Then the pressure difference driving the flow is the hydrostatic pres- 
sure outside the tube at y=h, less the hydrostatic pressure in the tube 
at this depth, 


Equation (1) gives an internal pressure difference of 


h h 
Pi(h)- p:(o) =\ dpi = ig + (Fo) sth) (38) Pli)/ dy (11) 


and an — — difference of 
Po(h) - Pe (0) = (fe & =) Sla)+ (3 ) 4} dy 


Since both pressures are atmospheric at the surface, subtraction of 


(12) 


equation (12) from equation (11) gives the pressure which drives the 


.-: {ffs s(h)] 36 + +Jeq)- -aly)3e f dy | (13) 


Dénsity is constant to within one per cent, and from the continuity 
equation, @v=constant. This implies that v*constant and, in fact, 


Will be assumed constant. Then the wall friction,7, is also constant: 


2 
= 3 £(R) Cae, (14) 
Here, £(R) is called the "resistance coefficient," (Rouse [5]). Due 


to the small value of viscosity for sea water, the flow will be as- 


sumed to be turbulent, and the relationship 


11 


4 
for smooth tubes will be used. Here, R is the Reynolds Nunber, 


R=2vr/Qy (16) 


where \ is the kinematic viscosity of sea water. 
Since the water in the tube is not accelerating, the driving 
force due to pressure must equal the retarding force due to wall 


friction: 
/\ Zee (17) 
ce @ 
Combining equations ie and (17) gives 
2 
Ve L F(R) = (18) 
ero 


ak 


Solving explicitly for Gp in equation (13) necessitates prescribing 








a salinity-depth curve. A fourth-power salinity profile 
S(y) = Q- id (19) 
gives the salinity curve with depth the same physically reasonable 
form as the temperature curve, The surface salinity is Q and 
P=(So-S,) pn’, ite (13) becomes 
g 2¢/( -ehah io, “Ap 
Ap = ghSt ee Sawi "-S(n]- (20) 
—— ih _ eAgh) 
de fexm i Gs (J-e4%h) - A, (I- e™ Ms 


430 A, A2 e ih. eAvzh 
Equation (20) is difficult to evaluate in its present form, but 


there are several approximations which can be made to simplify’ the eXe- 


pression without sacrificing accuracy. 


12 





First, assume that the rate of flow will be less than ten meters 
per second, Calculated values were much less than this value, and a 
velocity of ten meters per second would require an extremely large 
pressure force to drive the flow. 

Second, assume that the tube thickness is very much less than the 
inner radius of the tube. A thin tube is required to allow the water 
in the tube to reach approximately the same temperature as the water 
outside the tube as it flows upward. 

Third, assume terms that are 107° (order of magnitude) as large 
as the dominant tems are negligible. 

From equation (5), A,, is a function of < and @ and will always 
be a negative number, If h is greater than (say) 50 meters, (it will 
usually be much greater than this), A,h is a large negative nunber, 
Bide“! is noditmpie. 

The constants D, E, F, and G are such small numbers that, even 
when multiplied by 1, h, hé, and h? respectively, they are still negli- 
gible compared to other tems, 


Equation (20) becomes 


Ap = 9h3 SJa- = "=8(h)] (21) 


and combining tae (18) and (21) gives 


Phi 


This is a closed form of the solution for the flow velocity which in- 
Cludes the effects of vertical heat conduction. 
If vertical heat conduction is not considered, the governing 


differential equation (2) simplifies to 


3 


D (@-@) 
—t (eThe?dy ) =m £14 dy 5 (23) 


which has the solution 
h 


Oly) = yer fe*Tecr)a + ©lh) o~ tln-4) (24) 
4 
where 3 =2nL/wevrhe . 


Solving equation (24) for the fourth-power temperature profile gives 


0(y)= A-Cy'- y Cy? 12Cy? = oy 24C (84 +1) 
- 3 Z 8 (25) 
on thea) [#ch? | 12Ch? , 24C Gh )) 
b Se ¥ 
Using the fourth-power salinity profile and equation (25) for @(y), 


equation _ becomes 


Ne: We {Ce ph! -sts}- a/b + 2 24 ‘a 


where negligibly small terms have been ignored. “ solution for fourth- 
power temperature and salinity profiles when vertical heat conduction is 


not considered 1s 


y eh / \3e Ph” cp) 3? Ch’ 2H) 
v6.65 10 ela “4A Fe|&" F “430/% 7" ¥F o eB 
Equations (22) and (27) were solved using the temperature and sali- 
nity profiles pictured in figure 1, using identical tube dimensions, 
The flow velocity which was influenced by vertical heat conduction was less 
than one-tenth of one per cent greater than the flow velocity calculated 
ignoring vertical heat conduction (Table 1.). 
2.2 Solutions for linear temperature and salinity profiles. 
In order to consider the effect of different temperature and sali- 


nity profiles, linear profiles were assumed: 
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e(y)=M-Ny (28) 


and 


S(y) = Q@- Py. (29) 


The solution to equation (4) then becames 


A h_ |- earh 
O(¥) = «| fe me +] ise 1 M-Ny -X (30) 


where X=N Bu ® 
Solving equation (13) using the temperature profile of equation (30) 


gives 


9 Nigmen 
Ap= sre fa- 5th] + gh3Sx +456 * a me 


where negligibly small terms have been ignored. The solution for linear 
temperature and salinity profiles, including vertical heat conduction 


effects is 


PAIVANY, de}, _ Ph_ ee 
Hacsen* (3) |b a z ~ Sthif-gh3e x (32) 


af (A, -A1) - 


Assuming linear temperature and salinity profiles, equation (24) becomes 


<n fle -¥h-¥) i} + M- Ny m 


and the solution to equation (13) is 


Ap 2 ghifa- B- -s(h)}- 35 4-8] (34) 


where negligibly small terms have been ignored. Thus, the solution for 
flow velocity for linear temperature and salinity profiles where verti- 


cal heat conduction is ignored is 


IES 


Pees w/a pe fa. Pb sth] +93 ih MIAO) - 


Equations (32) and (35) were solved using the temperature and sali- 
nity profiles pictured in Figure 2, with identical tube dimensions, Again, 
the rate of flow influenced by vertical heat conduction was less than one- 
tenth of one percent greater than the rate of flow calculated ignoring 
vertical heat conduction. (Table 1.). 

It 1s interesting to compare the actual calculations of flow veloci- 
ty for the profiles in Figure 1 with those related to Figure 2. Vertical 
heat conduction was included, but this is not pertinent to the discussion 
below, The velocity solutions, equations (22), (25), (32) and (34) all 
contain the term [a-Ph/2-S(h} which is the average salinity of the water 
between the surface and the depth y=h, S, less the salinity at y=h, S(h). 
In making these calculations it was noted that the main factor causing 
the one centimeter per second difference in rate of flow between the pro- 
files in Figures 1 and 2 was the value of the S-S(h) tem, 

Calculating the value of S-S(h) for the salinity profiles in Figure 
1 and in Figure 2 showed that the two values were almost identical. 
However, the value of S obtained by applying the fourth-power salinity 
profile approximation to the salinity curve of Figure 1 gave an erroneous 
value for S-S(h). The error in the calculation of S by substituting an 
approximate function for the actual salinity curve caused a marked error 
in flow velocity result. This error was later verified by the computer 
program determination of flow velocity which required no approximation of 
‘S. This indicates the great importance of the salinity profile which 


will be discussed in more detail below, 
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The flow rates for the 8-S profiles contained in Figures 1 and 2 were 
also calculated utilizing a digital computer (see below). Here the ef- 
fects of vertical heat conduction were not included. The rate of flow in 
each caSe was within one per cent of the other, and was essentially the 
same as the flow rate calculated by hand from equation (32). 


TABLE 1 
RESULTS OF CALCULATIONS FOR RATE OF FLOW 


Case 8(y) & Sty) "Fixed # Hand ## Computer ## 
No, Profiles --Parameters" Calculations Calculations 
, Figure 1 Approximate 17.32 
2 Figure 2 Approximate 17,36 17,30 
3 Figure 1 Approximate lornsc* 

“ Figure 2 Ree amate 28.52 27,88 
5 Figure 1 Precise 17.40 
6 Figure 3 Precise 22,68 
7 Figure 4 Precise 26,60 
8 Figure 5 Precise 20.89 
9 Figure 6 Precise 8.03 
10 Figure 5 Precise*** 30.94 
11 Figure 3 Precise*** 33459 


* The flow velocity result is for the hand calculation of the fourth- 
power approximation of temperature and salinity profiles illustrated in 
Figure 1, 


** The inner radius of the tube was increased to 20 centimeters and the 
thickness of the tube increased to two centimeters. 


*k* The tube was inserted vertically so that the length of the tube equals 
the depth of the bottom of the tube. 


# Appendix E contains "fixed parameter" values, Except where otherwise 
noted, the inner radius of the tube is 0.1 meters and the wall thickness 
is 0.002 meters. 


## Flow velocity results are in centimeters per second. 
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Marked differences in rate of flow which occur when the tube is 
inserted into water which has different temperature and salinity pro- 
files are a direct result of the difference in the salinity profile fron 
one location to another. The temperature profiles did not significantly 
influence the rate of flow within the small but realistic range of 
temperature and salinity used by the author, 

2.3 Effects of vertical heat conduction, 

The importance of vertical heat conduction was investigated in 
another manner. The wall thickness of the tube was increased by a factor 
of ten. The flow velocity decreased, but the decrease was less than one 
per cent. Increasing the inner radius of the tube from ten centimeters 
to 20 centimeters and retaining the thicker-walled tube resulted in a 
Significant increase in rate of flow of 37%, 

Comparison of the results of the above calculations indicated that 
the vertical heat conduction due to vertical turbulent mixing was suf- 
ficient to warm the water in the tube. In order to determine the extent 
of the warming due to vertical heat conduction, the velocity of flow was 
then determined for the situation in which no horizontal heat conduction 
was allowed. In this case, the temperature in the tube is a function 
only of vertical heat conduction. The governing differential equation (2) 


Simplifies to 


D | Jo 
aE eT? dy = Kitr?dy D4 


(36) 





which has the solution 


O(y) = 20g e%- ©(0) - eee (37) 


whe re a=KL/cvhe ‘ 
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Solving equation (13) using equation (26) and linear temperature and 
Salinity profiles, 

4 delNh*, wh } 
vfoesu' ik “ -gh< Q- 2 st) “al ea =O 038) 
where negligibly small terms have been onitted. 

The resulting flow velocity is nearly identical to the flow velocity 
calculated using the same tube, but allowing horizontal heat conduction. 
The process of heat conduction in the vertical by turbulent mixing would 
require a longer period of initial pumping to reach a steady state con- 
dition, but the results indicate that the vertical heat conduction pro- 
cess is sufficient to heat the water inside the tube, 

The initial tube placement was at an angle of depression of 30 degrees 
to the water surface, This was done so that the water in the tube would 
require a longer time interval to reach the surface and undergo increased 
warming. Since the water in the tube is sufficiently warmed by vertical 
heat conduction to be less dense than the water surrounding it outside the 
tube, the tube placement was changed to a vertical position. The decrease 
in tube length resulted in a corresponding decrease in the effect of wall 
friction and a 48% increase in rate of flow. Such an increase is certainly 
Significant and will increase the volume rate of flow, 

3. Computer solutions. 

In order to determine the rate of flow for a "perpetual salt fountain" 
Situated in the ocean, actual oceanographic temperature and salinity data 
for a given area must be known. A program to obtain flow velocity so- 
lutions which neglected vertical heat conduction was written for a CDC 
1604 digital computer using FORTRAN 60 programming. Using this program, 
flow velocity was easily determined for any area in which temperature and 


Salinity data were available, 
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Temperature and salinity values for each one-meter interval over the 
entire length of the tube, tube dimensions, and ''fixed parameters" are en- 
tered on data cards in the program, "Fixed parameters" include kinematic 
viscosity of sea water (.)), density of sea water (OC), specific heat of 
sea water (c), acceleration of gravity (g), rate of change in density with 
change in salinity (°P/aS) and rate of change in density with change in 
temperature (PPAe). Appendix E contains the values used for "fixed 


parameters,"' 


Subroutine INTE computes the temperature profile in the tube, and 
subroutine DELTP computes the pressure force which drives the flow. Sub- 
routine INTE provided some interesting problems due to the size of the 
exponents of the exponential functions. A discussion of these problems, 
and the solution to them is included with the computer program in 
Appendix A. Subroutine DELTP is a trapezoidal integration which is as 
accurate as can be obtained from the small number of data points gene- 


rally given in a hydrographic sounding, 


Initial data used in the program was set up to treat the same case 
aS was calculated by Groves, Temperature and salinity profiles were 
devised so that values of @(0)-@(h), S-S(h), and the "fixed parameters" 
would be identical to those used by Groves. Additional test cases were 


run to compare with the results of hand calculations, 


Values for the constant tems %%e and*°/< were then calculated fron 
the equation of state given by Eckart [1] (see Appendix B) and more pre- 
cise values for the other "fixed parameters" were obtained and used with 
data from actual hydrographic soundings. Several cases of actual data 
were employed to get an idea of the range of values of flow velocity which 
could be expected in the ocean, The results are contained in Table 1, and 


the temperature-depth and salinity-depth curves are pictured in Figures 1-6, 
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In the actual ocean data cases, the depth to which the tube was ex- 
tended was determined by the depth of the halocline. The tube was always 
extended into the water just below a negative halocline,. 

4, Results. 
4.1 Vertical heat conduction versus no vertical heat conduction, 

Initial flow velocity calculations were made by hand using the initial 
values of the "fixed parameters" contained in Appendix E. The rate of 
flow for identical tube dimensions was calculated for @-S profiles con- 
tained in Figures 1 and 2. 

The calculation of flow rate was made twice for the linear @-S pro- 
files of Figure 2, once including the effects of vertical heat conduction 
and once neglecting the effects of vertical heat conduction. In each 
case, inclusion of vertical heat conduction had a negligible effect on 
the rate of flow. 

Flow velocity calculations were then made for the 8-S profiles of 
Figure 1 using the fourth-power temperature and salinity profiles. These 
functions were selected because they approximated the physical appearance 
of the curves in Figure 1 and because they were easily integrable. Flow 
rate solutions for the vertical heat conduction and no vertical heat con- 
duction equations were again identical to four significant figures. 

Flow velocity was then calculated for a situation in which no hori- 
zontal conduction of heat was allowed, so that only vertical heat con- 
duction was possible. The rate of flow for this calculation was nearly 
identical to the rate of flow calculated for two other cases which used 
the same tube dimensions, "fixed parameter" values, and 9-S profiles. 

The first allowed both horizontal and vertical heat conduction, and the 


second allowed only horizontal heat conduction. 
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On the basis of these calculations, it was concluded that the effect 
of vertical heat conduction on the rate of flow of the water inside the 
tube was negligible in the presence of horizontal conduction, Thus it 
was decided to make use of the computer program which did not include 
the effects of vertical heat conduction for the remaining calculations. 
4,2 Computer results versus Groves’ result. 

Utilizing the computer program allowed exact values of temperature 
and salinity at each depth to be used, Temperature and salinity values 
from Figures 1 and 2 were used in the program with a resulting rate of 
flow of 17.30 centimeters per second (cps) for Figure 1 profiles and 
17.32 cps for Figure 2 profiles. The 17.32 cps flow velocity associated 
with the profiles of Figure 2 was in close agreement with the rate of 
flow calculated by hand for the same profiles, The 17.30 cps rate of 
flow calculated for the profiles in Figure 1 was much greater than the 
rate of flow calculated by hand for the same profiles. Investigation of 
the hand calculation revealed that the S-S(h) term calculated by hand 
differed significantly from the actual value of S-S(h) obtained from 
Figure 1. This difference was attributed to use of the fourth-power 
salinity profile as an approximation of the actual salinity profile 
when making hand calculations, S-S(h) calculated from the curve of the 
fourth-power salinity profile is significantly different from the true 
value of S-S(h) in Figure 1, The influence of the approximation to the 
form of the temperature profile, however, was determined to be negligible. 

The profiles of temperature and salinity in Figures 1 and 2 were 
constructed so that the ani of S-S(h) and 9(0)-@(h) would be identical 
to those selected by Groves, Tube dimensions and "fixed parameter" 


values were also those used by Groves. Groves' calculated rate of flow 


dd 





was 17,50 cps as compared with a rate of flow of 17.32 cps obtained from 
the computer program. The difference in flow velocities can be attribu- 
ted to a rounding-off procedure used by Groves to simplify calculations, 

The values of she and 9e/y 5 calculated in Appendix B were then 
inserted into the original program along with more precise values of the 
other "fixed parameters" to determine what effects the more precise values 
would have on the rate of flow. The resulting velocity increase of eight 
centimeters per second was not significantly greater, but the more pre- 
cise values for the 'fixed parameters" were retained for the remaining 
temperature and salinity profiles. 

Temperature and salinity data from actual hydrographic soundings were 
also used, Figures 3 through 6 illustrate the actual profiles used, It 
was verified that areas of most pronounced negative haloclines produced the 
greatest rates of flow. 

A marked increase in rate of flow was noted when the tube was in- 
serted vertically instead of at a depression angle of 30 degrees to the 
ocean surface, The 48% increase in speed of flow will provide an equal 
increase in volume rate of flow. 

Rate of flow was also increased by increasing the inner radius of the 
tube. This is doubly significant since an increase in flow velocity in- 
creases volume rate of flow, and increased tube radius will also increase 
the volume rate of flow. 

9. Conclusions and acknowledgements. 

The test cases discussed in this paper have demonstrated some interest- 
ing features. First, the temperature-depth structure is of relatively 
little importance. The surface layer must be sufficiently warm to in- 


sure stability, but the magnitude of the temperature change between the 
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surface and the bottom of the tube has a negligible effect on the velo- 
city of flow. 

Second, the primary parameter influencing the rate of flow of water 
in the tube is the difference between the average salinity in the layer 
through which the tube extends, and the salinity of the water entering 
the bottom of the tube. To calculate the rate of flow for a given case 
to within one per cent of its precise value, a salinity-depth curve must 
be known. A simple area-averaging process can be used to detemine S, 


the average salinity, and S(h) is known, Then, 


yy, be 
vifpesnot Sa] = g hho S-str]e ow 


Third, equation (39) gives excellent results because the effects 
of vertical heat conduction heat the water in the tube sufficiently to 
make it buoyant. Although vertical heat conduction does not directly 
increase the velocity of flow, the warming effects allow the substi- 
tution of a poorer heat-conducting material for the copper tube. A 
metal which is resistant to sea water corrosion would be a better material 
for the tube. In addition, a relatively thick tube with a larger inner 
radius may be used, This will increase the volume rate of flow and 
provide a more durable tube. 

For example, using data from the equatorial Atlantic Ocean region 
in case #11 (Figure 3), the volume rate of flow is 10.5 liters per second. 
This is a 90% increase in volume rate of flow over that calculated by 
Groves. When the inner radius of the tube was enlarged in case #4 
(Figure 2), the volume rate of flow increased to 34,6 liters per second. 
This is an increase of 530%, which is definitely significant, 

It is therefore concluded that the "perpetual salt fountain" is more 


than an “oceanographic curiosity."' Upward flow of nutrient-rich water 
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from well below the surface in sufficient vclume could make a fertile 
fishing area out of an area which is, at present, barren. A relatively 
inexpensive and simple method for providing this nutrient-rich water, 
such as the method described in this paper, is worthy of additional 
experimentation. 

The author gratefully acknowledges the guidance, encouraging 
assistance and critical comments of Assistant Professor T. Green. My 
thanks also go to Mr. R. D. Brumell of the United States Naval Post- 
graduate School's canputer facility staff for his valuable programming 


assistance, 
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APPENDIX A 


COMPUTER PROGRAM FOR FLOW VELOCITY CALCULATIONS 


ee JOBO532FMsHINMANsKeGe s9FLOWUP 


5) 0 


10 


eZ 


i2 


2D 
36 
20 
At 
22 
st 


57 


40 
41 
43 


44 
45 


42 


PROGRAM FLOWUP 

DIMENSION 1T(401)5 S(401).s PHI(401) 
COMMON NVT »HsGAMMA sDELHsPHI 9S5TsSALTHs THETHs Gy DELTP » CHDS » CHDT 
READ 50,5 NVT 

FORMAT (1215) 

READ 109 (T(I)sS(I)sT=l1sNVT) 

FORMAT (6E12¢6) 

READ 109 ReExsCONsPLoWsDsSHsDEPTHsV0sEP1l»sTHETH 
READ 109 HsGsDELHsSALTHsCHDS »CHDT 
PRINT 129 ReEsCONoPLoWoD 

PRINT 12s SHsDEPTH»sVOsEP1sTHETH 
PRINT 129 HsGsDELHsSALTHsCHDS ,CHDT 
FORMAT (6E12-26) 

V=V 

V=V+0.01 

RE=2eO*VE¥ER/SE 

GAMMA=4 eO*CON¥PL/ (W¥D¥SH*¥EXREXDEPTH) 
CALL INTE 

CALL FINAL 

IF(DELTP-0e0) 555955956 

V=V-0.01 

GO TO 57 

TAU=DELTP#¥R/(22e0*PL) 

IF(RE-200020) 20920921 

FUNC=64e0/RE 

GO TO 22 

FUNC=023164/( (RE) ¥*¥#( 120/40) ) 
VEL=(8eO0*# TAU/(FUNC#D) )**¥(120/2.0) 
FORMAT (30H THE UPWARD FLOW VELOCITY IS 2E20. 6) 
IF(VEL-V) 55924919 

V=V+0.001 

RE=2eO0*V¥R/E 

GAMMA=4 eO*CON¥PL/ (W#D¥SH¥E#REXDEPTH) 
CAL Ear e 

CALL FINAL 

IF(DELTP-0e0) 40940941 

V=V-0.001 

GO TO 42 

TAU=DELTP#R/(2e0*PL) 

IF(RE-200020) 43943544 

FUNC=64.0/RE 

GO TO 45 

FUNC =02¢3164/((RE)**¥(0025)) 
VEL=(8e0*TAU/(FUNC*D) )**(025) 
IF(VEL=V) 40924557 

V=V+0.0001 
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10 


RE=2-0#V#¥R/SE 

GAMMA=4 eO*CON*#PL/ (W#D*¥SH*E*#REXDEPTH) 
CALL SINE 

CALL FINAL 

IF (DELTP~0-0) 46946947 

V=V-0e0001 

GO TO 7¢ 

TAU=DELFP#R/(220*PL) 

IF(RE=2Q900-0) 71ls71ls72 
FUNC=640e0/RE 

GO T0723 

FUNC=003164/ ( (RE) ¥¥(140/400),) | 
VEL=(8e0*TAU/ (FUNC*D) ) ¥¥( 120/220) 
IF(VEL=$V) 46324942 oe 
V=V+0.00001 ns 

RE=2eO*V*¥R/E 

PRINT 60» RE | 

FORMAT (25H THE REYNOLDS- NUMBER IS E2026) 
GAMMA=4 eO*CON#PL/ (W¥D*¥SH*® EXREXDEPTH) 
PRINT 619 GAMMA 

FORMAT (14H GAMMA EQUALS E20¢6) 
CAEL INTE 

CALL FINAL 

PRINT 63s DELTP 

FORMAT (14H DELTP EQUALS E20e6) 
TAU=DELTP¥R/(2-0*PL) 

PRINT 629 TAU 

FORMAT (12H TAU EQUALS E20-6) 
IF(RE=2000¢0) 745974575 

FUNC=64 e-O/RE 

GO TO 76 
FUNC=023164/( (RE) **(1¢0/4e0) ) 
VEL=(8eO*XTAU/ (FUNC#D) ) ¥¥(12¢0/2.0) 
PRINT lls VEL»V 
DIFF=ABSF((VEL-V)/V) 

IF(DIFF+EP1) 24924570 

STOP 

END 

SUBROUTINE INTE 

DIMENSION T(401)9 S(401)s PHI(401) 
COMMON NVT»H»GAMMA sDELH »PHI 9ST sSALTH»s THETHsGsDELTP » CHDS» CHDT 
DO 10 J=l»sNVT 

PHI (J)=020 

CONTINUE 

ITJ=NVT 

PHI (IJ)=THETH 


. BIG=GAMMA+1.0 
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BIGEX=E XPF (GAMMA) 

BNE =GAMMA=1.0 

ITJ=IJ-1 

K=ITJ+1 
XI=(TCIS)—-(T(K)*BIG) 7 (GAMMA*BIGEX) 


i 


o10, 


40 


DO 


XJ=((TC IS) *¥BNE)+T(K) ) “GAMMA 

XK=PHI(K) /BIGEX 

PHI (IJ) =XI+XJ+XK 

IF (IJ-1) 30530220 

RETURN - 

END } 

SUBROUTINE FINAL 

DIMENSION T(401)»9 S(401)s PHI(401) 

COMMON NVT»sH»GAMMA sDELHsPHI 9Ss9T»SALTHs THETHsGsDELTP »CHDS » CHDT 
TOTAL=0 0 
[=] | , 
PRESS=(S(1)—-SALTH) *CHDS+(T(1)—PHI(1))*CHDT ; 
l=1+] | 

SUMP=((S(1)-—SALTH) *CHDS)+((T(1)—PHI(1)) *CHDT) 
TOM=DELH* (PRESS+02e5*(SUMP-PRESS)) 

PRESS=SUMP 

TOTAL=TOTAL+TOM 

IF (NVT-I) 50950940 

DELTP=G*¥TOTAL 

RETURN 

END 

END 


28 








The computer program is an iteration process for obtaining the 
velocity of the water in the tube. It was written to solve actual 
oceanographic cases using the form of the solution given by Groves [2]. 


Subroutine DELTP is utilized to compute the pressure integral, 


h 
[p= af {[su- stn] Se +[eq)- oly)) 38 f dy ee: 


The subroutine is a trapezoidal integration which assumes a linear 
relationship between successive points. No greater accuracy is desirable 
because hydrographic soundings contain salinity and temperature data for 
relatively few depths, and the linear relationship is generally assumed 
when plotting the original data, 

Subroutine INTE is designed to canpute the temperature of the water 


at each one-meter interval of depth inside the tube using 


O(y) = se*4] OG)e "dT + ohbe” cha) (A2) 
¢ 


Equation (AZ) relates water temperature in the tube to depth, but it 
cannot be solved in this form on the computer since the capacity of the 
CDC 1604 is exceeded by exponential functions with exponents of magni- 


tude greater than 709. A different form of equation (A2) was necessary 


for the computer program. h 
Let z(y) = ge%| ©(T) oe eaT (A3) 
md ¥(y) = O(n) @ ~Mh¥) | a 


mm Oa) = 24) + YC4)> a 


In order to get Z(y-Ah), we let 


X (y-Ah) = ye" an) foene 


Zo 


"dT (A6) 


and then from equations (A3) and (A6) 

z(y-ah)= X(y-ah) +e ZG). 
We now assume that 0(T), where T is a "dummy" variable of integration, 
varies linearly between consecutive one-meter intervals, and limit the 
interval, Ah, to one meter. 

Then 


zi -Ah) - te *ls1eG)- e(y-ah)} -|el +(4-1)(y- if 


Canbining iar (AS), (A7), and (A8) 


PCy -Ah) = i “tely)G+!) 7" 6(4-4h)| 
ste [ ely) + (4-1) e(y-ah)] § * edly). (A9) 


Since $(h)=0(h) is a boundary condition, the first %(y- Ah)will be cal- 


culated for a depth of (h-1) meters, This integration subroutine will 


function so long as §< 709. 
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APPENDIX B 
COMPUTATIONS OF 4g ann °S% 5 
The values for the rate of change of density with changing tempera- 
ture and salinity were assumed constant in the solution of the pressure 
integral, Additional investigation of the value of these constants was 
considered important. 


Assuming that the effects of pressure are negligible, 


oe d 
e(S,@) = Ri + (s-s,)($$), r(e-e)(SS) + on a (B1) 
Following Eckart [1], 


(p+ P)(A.-A)=K (82) 


where 

A = jot (i- t) (B3) 
and 
ck AO ea +0.375T% +38 (B4) 


K = [17 195 +0, N25 T-O,000745 74 - (0,0380 +0,00017) S ]x 0" 95) 


may. = Cons tain. = 5a2@). (B6) 


Conbining equations (B2) and (B3) gives 


e= (y -[101 (Ag — tat. (B7) 


Taking the partial derivatives of Q with respect to 9 and S in 
equation (B7), and evaluating the result at an average temperature, 
Q=15°¢C, an average salinity, S=34°/o0, and a zero pressure, gives 


I°f S =0. 77667kg/m>* o/ooand Ko =(), 19987kg /m>* Oc. 
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APPENDIX C 


WALL FRICTION CALCULATION 


It was considered that the effects of the stress on the water by 
the wall of the tube might cause significant heating of the water, The 
work done on the water per unit mass per time is 

T 20 4 dy Vv DD ae | 

————_ = —— joules 

eT t*dy e ( g*Sec ), (C1) 
The work per unit mass is 


La Ly a ar ly joules (C2) 
er Vere ( 8 ) 


Dividing by the specific heat, c, and converting joules to calories using 








0.2389 calories/joule at 15°C, 
ct. (0.2389) = AT (*). (C3) 
eC A 
Applying the data used in the solution of case #1, the temperature change 
of the water in the tube due to wall friction is 0.00025°C for a velocity 
of 17.32 cps and a pipe length of 600 meters, This change in temperature 


is negligible. 
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APPENDIX D 
COEFFICIENT OF VERTICAL HEAT 
CONDUCTION FOR TURBULENT FLOW 


The value for the coefficient of vertical heat conduction for turbulent 
flow, K, was investigated, From Schlichting, [6], the coefficient of eddy 
conductivity of heat, Ag» is at least as large as the coefficient of eddy 
viscosity, Ay, over the entire pipe. We assume the turbulence to be ap- 
proximately isotropic, so that turbulent transfer along the pipe is the 
same order of magnitude as that in any cross section, 

A good estimate of the effective value of Aq can then be obtained by 
calculating an average of Ayover a cross section. This is easily done 


uSing mixing length theory. The result is 


ar 
- > e 2 im Oe 
Ay = Ar = S| Uy# am(r-2z) dz — re (D1) 


where r is the interior radius of the tube, z is measured inward from the 
wall, ux is the "shear velocity" (see Schlichting), and k=0.4, 
Then, 


K = Ag Cp, (D2) 


This has a typical value of 5.0 cal./m*°C*sec. for the rates of flow 


and tube size considered, 
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APPENDIX E 


"FIXED PARAMETER" VALUES 


The expression ''fixed parameters" refers to values for kinematic 
viscosity of sea water (\\), conductivity of tube material (n), densi- 
ty of sea water (QO), specific heat of sea water ( c ), acceleration of 
gravity ( g ), rate of change of density with changing salinity (20/35) 
and rate of change of density with changing temperature (°P/e). Two 
sets of values for the ''fixed parameters" have been used. One Set con- 
tains rounded-off values for hand calculations, and the other set contains 
more precise values for the computer program, One could then determine the . 


effect on rate of flow of changes in the values of the ''fixed parameters," 


APPROXIMATE PRECISE UNITS 
) 107° 1,22x1079 m2/sec 
rr (copper) 10° 10° cal/m*°c,sec 
Q 10° 1.026X10° kgm/m> 
C 10° 9, 33X10" cal/kgm+°C 
yY 10 9.8 m/ sec’ 
90/35 0.7 0.77667 kgm /m? 
-0,2 -0.19987 kgm /m>* °C 


de/JO 
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